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摘 要 本 文 基于 多 源 X 射线 层 析 成 像 技 术 及 联合 代数 重建 算法 (SART) 首次 重 构 出 黏 性 B 类 颗粒 流 化 时 局 部 气泡 的 
三 维 形态 , 并 结合 压力 脉动 分 析 技 术 获 得 了 丰富 的 全 局 流 态 化 信息 . 实验 结果 表明 ,， 黏 性 力 的 增加 强化 了 气泡 间 的 合并 作 
用 , 导致 气泡 的 尺寸 上 升 , 频率 和 速度 下 降 ， 当 床 温 超过 35°C 后 ， 节 消 状 态 被 激发 ， 并 与 自由 流 化 状态 交替 产生 . 而 且 ， 
黏 性 力 越 大 , 节 涌 在 流 化 过 程 中 所 占 的 时 间 比 例 越 长 . 


关键 词 竺 性 颗粒 ; X 射线 层 析 成 像 压力 脉动 分 析 ; 流 态 化 
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大 量 颗粒 的 能 力 ， 在 不 同 工 业 领 域 中 得 到 了 广泛 能 
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Abstract This paper first reconstructed three dimensional version of local bubbles in a fluidized bed 
processing cohesive Geldart B particles with a combination of X-ray tomography and Simultaneous 
Algebraic _ Reconstruction Technique (SART). Additionally, global fluidization dynamics were also 
obtained using in-bed pressure fluctuation analysis. The experimental results show that the preserve 
of cohesive force facilitates bubble coalescence, thus increasing the bubble size and decreasing the 
bubble frequency and bubble rising velocity. When the bed temperature exceeds 35°C, slugging takes 
place alternatively with normal fluidization. Moreover, the duration of slugging sharply increases with 
the cohesive force. 


Key words cohesive particles; X-ray tomography; pressure fluctuation analysis; fluidization 


流 化 床 反应 器 因 其 高 效 的 传 热 传 质 与 连续 处 理 ” 流 态 化 特性 研究 起 步 较 早 ， 最 初 集中 在 湿 颗 粒 的 临 
为。 界 流 化 特性 上 3。 近年 来 ， 研 究 重点 逐渐 向 流 型 转 


立 用 


. 伴随 着 流 化 床 处 理 对 象 与 手段 的 多 样 化 , 颗粒 ”化 67、 乳化 相 流动 形态 加 及 流 化 失效 判定 四 等 方 


间 经 常 出 现 黏 性 力 (如 干燥 、 造 粒 、 涂 层 等 过 程 ), 导 i 迁移 . 需要 指出 的 是 , 从 实验 角度 出 发 , 在 颗粒 间 
致 颗粒 的 流动 性 降低 ， 流 态 化 行为 发 生 改 变 ， 从 而 ”均匀 引入 液 桥 力 的 难度 较 大 ， 导 致 条 性 力 易 呈 各 向 
影响 流 化 床 内 
因此 ， 揭 示 符 性 力 对 流 态 化 特性 的 影响 规律 和 内 在 ”等 采用 “多 聚 物 涂 层 " 法 引入 固 桥 力作 为 接触 式 黏 性 


的 反应 效率 以 及 反应 器 的 稳定 运行 。 异性， 降低 结果 的 准确 性 . 针对 此 不 足 ，Shabanian 


证 到 


机 理 ， 对 流 化 床 反 应 器 的 设计 和 运行 具有 重要 的 指 。 力 bo, 并 系统 地 研究 了 不 同 务 性 力作 用 下 流 化 状态 
导 意 义 ， 从 鼓 泡 流 化 向 清流 流 化 转变 特性 的 差异 552 

黏 性 力 的 种 类 众多 ， 本 文 着 重 探讨 接触 式 攻 性 由 上 可 知 ， 已 有 研究 多 集中 在 黏 性 力 对 乳化 相 
力 ， 即 在 颗粒 接触 后 产生 的 吸引 力 中 ,其 中 以 液 桥 。 流动 行为 的 影响 规律 上 , 对 于 气泡 相 的 研究 较 少 , 而 
力 外 和 固 桥 力 加 最 具 代 表 性 . 对 于 液 桥 力 存在 下 的 ”气泡 特性 恰恰 是 影响 鼓 泡 流 化 床 运行 中 传 热 传 质 的 


收 稿 日 期 ，2017-01-18; 修订 日 期 ，2017-07-23 


基金 项 目 : 国家 重点 研发 项 目 (No.2016YFB0600802); 中 央 高 校 基本 科研 业务 费 专项 资金 资助 (No.222201717004) 


作者 简介 : 马 吉 亮 (1984- )， 男 , 讲师 , 主要 从 事 气 固 两 相 流 动 特性 研究 , Email: jlma@seu.edu.cn 


chinaXiv:201710.00074Vv1 


ChinaXiv 合 作 期 刊 
8 期 马 吉 亮 等 : 笑 性 Geldart B 类 颗粒 流 化 特性 实验 研究 1703 
重要 参数 , 此 外 ， 流 态 化 在 高 共性 力作 用 下 的 失效 测 点 构成 . 为 了 避免 X 射线 辐射 对 实验 人 员 健 康 的 
机 理 ， 尤 其 是 失效 前 后 动力 学 参数 的 变化 规律 尚未 ”影响 , 所 有 实验 设备 均 布 置 在 由 5 mm 铝板 围 成 的 
见报 道 . 房间 中 . 实验 时 通过 电脑 和 数据 采集 卡 实现 对 系统 
本 文采 用 实验 手段 研究 黏 性 B 类 颗粒 在 三 维 鼓 。 启 停 和 实验 数据 采集 的 控制 . X 射线 强度 数据 的 采 
泡 流 化 床 中 的 流 化 特性 : 利用 “多 聚 物 涂 层 ”法 引入 集 频率 为 2500 Hz, 采集 时 长 1 min; 压力 脉动 信号 
接触 式 猪 性 力 ， 并 通过 调节 床 温 实现 对 黏 性 力 的 精确 。 的 采集 频率 为 1250 Hz, 采集 时 长 10 min. 
深 制 ; 基于 荷兰 代 尔 夫 特 理工 大 学 搭建 的 多 源 XX 射线 。 1.2 实验 材料 
层 析 成 像 系 统 重 构 床 内 的 局 部 三 维 气泡 形态 , 同时 结 本 实验 采用 “多 聚 物 涂 层 ” 法 引入 素性 力 9， 
合 压 力 脉动 分 析 技 术 提取 全 局 流 态 化 信息 ， 重 点 关注 ” 即 在 惰性 颗粒 表面 均匀 包 衰 一 层 多 聚 物 ， 利 用 多 聚 
流 化 失效 前 后 流体 动力 学 特征 的 变化 规 得 


物 茜 度 对 温度 的 高 敏感 性 在 颗粒 间 引 发 固 桥 力 即 接 

] 实验 到 钥 式 铬 性 力 。 该 方法 的 优点 是 重复 性 好 ， 医 性 力 分 
实验 系统 及 方法 布 均匀 ， 操 作 温 度 低 . 实验 采用 直径 0.6 mm, 密度 
和 2500 kg/m? 的 臻 璃 珠 作为 “基底 颗粒 "， 其 临界 流 


图 1 所 示 为 三 维 流 化 床 实验 系统 图 . 流 化 床 本 ”化 风速 Das 为 0.325 m/s， 属 于 Geldart B 类 颗粒 ; 
体内 径 140 mm， 高 1800 mm， 采 用 金属 烧结 板 布 ”采用 聚 乙 基 丙 烯 酸 酯 (PEA) 与 聚 甲 基 丙烯 酸 甲 酯 
风 . 实验 通过 电 加 热 器 及 PID 控制 器 调节 流 化 风 的 PMMA) 的 混合 物 作为 多 聚 物 涂 层 . 涂 层 过程 在 滚 
温度 ， 从 而 达到 控制 床 层 温度 及 颗粒 间 黏 性 力 的 目 。 简 式 包 衣 机 内 完成 ， 操 作 步 又 已 在 文献 [10] 中 详细 
的 。 床 层 温度 由 布置 在 距 布 风 板 110 mm 处 的 热电 ” 介绍 . 图 2 为 涂 层 前 后 颗粒 表面 的 SEM 图 . 本 文中 ， 
偶 进行 监测 。 在 距离 布 风 板 30 mm、110 mm、190 多 聚 物 涂 层 厚度 约 为 10 hm。 

mm 和 270 mm 处 分 别 布置 了 Kistler 7261 型 高 敏 
压力 传感器 对 床 内 表 压 进行 测量 . 同时 , 基于 X 射 
线 层 析 成 像 技 术 检 测 横 截 面 上 的 空隙 分 布 ， 其 原理 
是 当 X 射线 穿 过 颗粒 层 时 ， 射 线 强 度 会 发 生 指数 衰 
减 : T= Joe-*, 式 中 jo 和 了 分 别 为 穿 过 颗粒 层 前 
后 X 射线 的 强度 ,x 为 颗粒 层 长 度 ,/ 为 衰减 系数 . 
实验 时 ,通过 测量 每 条 X 射线 强度 的 衰减 率 即 可 获 


一 、 


得 该 条 射线 经 过 的 颗粒 层 长 度 .本 实验 采用 三 个 X 图 2 颗粒 涂 层 前 后 的 表面 特性 
射线 源 互 呈 120° 角 布 置 在 流 化 床 周围 ; 每 个 射线 源 (a) 涂 层 前 ，(b) 涂 层 后 
对 应 一 组 探测 器 ， 布置 在 流 化 床 外 | 与 身 线 源 相对 Fig. 2 Surface morphology of particles before and after 


coating (a) coated (b) uncoated 


的 位 置 ; 每 组 探测 器 分 为 上 下 两 排 , 每 排 由 32 个 探 


1.3 数据 处 理 方法 


六 本 文采 用 的 压力 脉动 分 析 算法 主要 包括 功率 谱 
bea pt 、 密度 和 平均 频率 等 常用 方法 ， 因 此 不 再 歼 述 . 下 面 
和 将 对 气泡 形态 重 构 算法 进行 简要 介绍 ， 

名 Wn 气泡 重 构 过 程 分 为 三 步 : 1) 将 3 组 探测 器 共 192 


个 探测 点 检测 到 的 X 射线 强度 代入 标定 函数 ,转化 

为 每 条 X 射线 经 过 的 颗粒 层 长 度 ，P; 2) 将 重 构 区 
上 层 探测 面 域 ( 流 化 床 内 壁 的 外 切 四 边 形 ) 网 格 化 ， 并 计算 射线 
在 网 格 内 的 权重 因子 W， 本 文中 重 构 区 域 的 像素 分 
辩 率 是 80x80; 3) 求解 式 (1), 获得 每 个 像素 内 的 颗 
下 层 探测 面 粒 份额 a. 


热电 偶 


N 
pi= >》 Wijay (1) 
j=1 


到 1 三 维 鼓 泡 流 化 床 实验 系统 式 中 , zi 是 射线 i 经 过 的 颗粒 层 长 度 ; Wi; 是 射线 i 
Fig. 1 Schematic of experimental setup 在 像素 用 内 的 权重 因子 即 射线 长 度 ; N 是 图 像 内 的 
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像素 数 . 2 结果 与 讨论 

本 文采 用 Anderson 和 Kak 于 1984 年 提出 的 联 
合 代数 重建 算法 (SART) 求解 式 (1)。 算法 的 实现 过 图 5 是 不 同 笑 性 力作 用 下 , 距 布 风 板 70 mm 横 
程 已 在 文献 [13] 中 进行 了 详细 的 论述 . 本 文中 , 重 构 ”截面 上 的 三 维 气泡 重 构图 . 当 床 温 为 30"C 时 ,经 过 
所 用 数据 为 每 10 组 原始 数据 的 平均 值 , 意味 着 数据 ” 横 截 面 的 气泡 体积 较 小 、 数 量 众多 旦 存在 彼此 相连 
的 有 效 采样 频率 由 2500 帧 / 秒 降 为 250 帧 / 秒 . 这 昌 的 多 重 气泡 结构 ; 随 着 黏 性 力 的 上 升 ， 气 泡 频率 逐 
然 降低 了 重 构 结果 的 时 间 分 辨 率 , 但 却 有 效 地 减 小 渐 降 低 ， 气 泡 尺 寸 逐渐 增加 ， 意 味 着 气泡 间 的 合并 
了 由 射线 源 强度 波动 造成 的 误差 . 图 3 对 比 了 上 下 ”作用 增强 ; 当 床 温 达 到 45°C 时 , 床 内 出 现 了 横向 尺 
两 层 探测 面 在 同一 时 刻 的 重 构图 像 ， 图 中 白色 区 域 。 寸 与 床 内 径 相 同 的 大 气泡 , 流 化 状态 转变 为 “ 节 涌 ”。 
被 气泡 占据 . 将 不 同时 刻 的 重 构图 像 进 行 琶 加 , 即 可 ”气泡 尺寸 与 黏 性 力 间 的 正比 关系 证 明了 随 着 黏 性 力 
获得 经 过 探测 面 的 增加 ,气泡 相 与 乳化 相 之 间 的 传 质 逐 渐 恶 化 。 


中 红色 对 象 为 气泡 。 
上 下 两 层 探测 面 上 出 现 的 时 间 差 估算 该 气泡 的 上 升 


的 三 维 气 泡 图 像 ， 如 图 4 所 示 ， 图 
包 . 时 


构 后 ， 可 以 基于 同一 气泡 在 


速度 . 同 理 , 也 可 以 针对 三 维 气泡 提取 尺寸 、 形 状 、 
频率 、 分 布 等 特征 参数 ， 


(b) 


图 3 上 下 两 层 X 射线 平面 重 构 结果 对 比 
(a) 上 层 探测 面 ,，(b) 下 层 探测 面 


Fig. 3 Reconstructed results for upper and lower planes 
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Fig. 
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图 4 上 下 两 层 探测 器 3D 气泡 重 构 对 比 
(a) 上 层 探测 面 ; (b) 下 层 探测 面 
4 Reconstructed bubbles for the upper and lower planes 


(a) upper plane (b) lower plane 


1490 ys 70140 
y/mm x/mm 
图 5 距 布 风 板 70 mm 处 的 三 维 重 构 气 泡 形态 随 黏 性 力 的 变化 

规律 , 静止 床 高 : 300 mm, 流 化 风速 : 2.5Cmf 


Fig. 5 Effects of cohesive force on the bubbles reconstructed 


14908e— 140 140 140 
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at 70 mm from gas distributor. Static bed height: 300 mm, 
fluidization gas velocity: 2.5Umf 


图 6 所 示 为 不 同 黏 性 力作 用 下 各 高 度 横 截 面 上 
的 气泡 分 布 热 图 . 图 中 各 像素 的 颜色 表征 横 截 面 上 
每 一 点 被 气泡 占据 的 概率 . 在 靠近 布 风 板 的 区 域 (图 
6a, b)， 随 着 黏 性 力 的 增加 ,气泡 的 合并 增强 ,， 尺 才 
增 大 ， 因 此 床 截 面 中 心 区域 的 气泡 高 概率 区 域 逐 渐 
六 张 . 随 着 测量 高 度 的 提升 ， 气 泡 概 率 的 变化 特 怕 
与 句 性 力 紧密 相关 : 当 和 素性 力 较 低 时 (30"C)， 气 泡 
的 横向 尺寸 伴随 着 气泡 合并 逐渐 增加 ， 气 泡 高 概率 
区 的 面积 增 大 ; 当 黏 性 力 较 大 时 (45"C)， 节 涌 产 生 
的 料 栓 完 全 占据 了 横 截 面 ， 导 致 气泡 出 现 的 概率 降 
氏 ,， 因 此 测量 高 度 越 大 ， 高 概率 区 的 面积 越 小 . 
图 7 为 不 同 高 度 截面 上 的 平均 气泡 上 升 速度 随 
黏 性 力 的 变化 规律 . 随 着 攻 性 力 的 增加 ， 平 均 气 泡 
上 升 速度 逐渐 降低 . 这 与 传统 的 气泡 理论 即 自由 气 
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泡 的 上 升 速度 与 尺寸 呈正 比 关系 [3 相 矛 盾 ， 意 味 
下 ， 床 内 存在 许多 大 尺寸 、 低 速度 
泡 ,基于 图 5 的 观测 结果 可 以 推测 ， 这 些 低速 
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气泡 是 伴随 着 节 涌 产生 的 气 栓 . 理论 上 ， 气 栓 的 上 
升 过 程 与 传统 气泡 不 同 : 由 于 气 栓 同时 受到 颗粒 重 


力 以 及 料 栓 与 边 壁 之 间 摩 擦 力 的 影响 ， 导 至 
速度 较 小 b4， 


符 随 着 屠 性 力 的 增 


其 上 升 


加 ， 节 涌 发 生 的 频 


率 提高 ， 气 栓 在 探测 到 的 气泡 中 所 占 比 例 增加 ， 导 
致 平均 气泡 上 升 速度 减 小 .需要 注意 的 是 ， 平 均 气 


泡 速度 在 床 温 超过 3 


是 床 内 节 涌 频繁 发 生 的 起 始点 ， 


5°C 后 大 幅 下 降 , 意味 着 该 温度 


和 A 30°C 40°C A 
3 
0.4 
证 | (b) 
叫 
尺 0.2 
电 ) 总 等 
(a) 
0.0 
> 
黏 性 力 
图 6 不 同 高 度 上 的 气泡 分 布 热 图 随处 性 力 的 变化 规律 


(a) 70 mm; (b) 150 mm; (c) 230 mm 
静止 床 高 : 300 mm, 流 化 风速 : 2.5Dmf 


Fig. 6 Contour plots of bubble distribution as a function of 


cohesive force at different heights. Static bed height: 300 mm, 
fluidization gas velocity: 2.5Umf. (a) 70 mm; (b) 150 


mm; (c) 230 mm 
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图 7 不 同 高 度 上 的 平均 气泡 上 升 速 


床 温 /*C 


止 床 高 : 300 mm, 流 化 风速 : 2.5Umf 


度 随 竺 性 力 的 变化 规律 ， 静 


Effects of cohesive force on the averaged rising 


velocity of bubbles at different heights. Static bed height: 300 


mm, fluidization gas velocity: 2.5Umr 


度 (PSD) 随和 共性 力 的 变化 规律 . 
寺 ，PSD 强度 总 体 呈 单 峰 分 布 ， 主 
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图 8 为 不 同 高 度 处 的 压力 脉动 信号 功率 谱 密 


当 床 温 为 30°C 


主峰 外 ，PSD 
峰 , 对 应 大 气泡 的 破裂 或 多 有 
性 力 较 大 


频 在 1 Hz 左右 。 


1 线 中 还 存在 若干 能 基 较 小 的 独立 
气泡 结构 的 产生 . 当 黏 
和 ,PSD 强度 由 单 峰 分 布 转变 为 双 峰 分 布 ， 


节 涌 持续 时 间 逐 渐 延 长 ， 


图 8 


功率 谱 密度 /kPa2/Hz 


其 中 低频 峰 出 现在 0.1 Hz 左右 . 由 图 5 可 知 ， 双 峰 
分 别 对 应 节 涌 和 正常 流 化 . 而且 由 于 正常 流 化 状态 
在 持续 时 间 上 占据 主导 地 位 ， 因 此 节 涌 对 应 的 频率 
仅 有 0.1 Hz 左右 . 随 着 黏 性 力 的 增加 ， 床 内 产生 了 
大 量 的 复合 结构 气泡 [加 ， 因 此 PSD 曲线 的 频带 宽 
度 逐 渐 增 加 ; 同时 ， 曲 线 整体 开始 向 低频 方向 移动 ， 
说 明 床 内 气泡 数量 逐渐 降低 ， 其 中 ， 低 频 峰 的 峰值 


随 务 性 力 的 上 升 而 增加 ,说明 气 栓 的 尺寸 逐渐 增 大 ， 


和 | 低频 峰 


2.0 


in 


= 


So 


全 
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-一 407C 


”频率 /Hz 


不 同 高 度 处 压力 脉动 信号 功率 谱 密 度 (PSD) 随 柯 性 力 的 
变化 规律 , 流 化 风速 ; 2.5Uut， 静止 床 高 : 300 mm 


Fig. 8 Power spectral density of pressure signal measured at 


different heights. Fluidization gas velocity: 2.5Umt, static bed 


height: 300 mm 


图 9 为 功率 谱 密度 的 平均 频率 f。 随 苛 性 力 和 
测量 高 度 的 变化 规律 。f。 定义 为 压力 信号 功率 谱 密 


图 9 所 示 , fs 随 


压力 信号 中 处 于 低频 区 的 “气泡 通 


说 明 此 效应 


度 图 中 ， 将 信号 能 量 二 等 分 的 频率 ， 反 映 了 不 同 物 
理 现象 在 整体 信号 中 所 占 比 例 。 如 

着 共性 力 的 增加 逐渐 降低 : 在 35°C 之 前 ， 气 泡 的 
合并 作用 增强 ， 压 

道 效应 ”在 信号 中 占据 的 比例 提高 中 ， 导 致 f。 逐 
渐 降 低 ， 但 由 于 fo 的 下 降幅 度 较 小 ， 

对 压力 信号 的 能 量 分 布 影 响 有 限 ; 当 床 温 超 过 35°C 


后 , 床 内 出 现 低频 节 涌 ， 
动 ， 导致 fo 继续 下 降 . 


味 着 
上 升 


促使 能 量 分 布 向 低频 方向 移 


而 且 ， 下 降幅 度 迅 速 增加 意 


节 涌 在 流 化 过 程 中 所 占 比例 在 35°C 之 后 迅速 
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图 9 香 性 力 及 测量 高 度 对 功率 谱 密 度 平均 频率 的 影响 , 流 化 风 


速 : 2.5Umf， 静 止 床 高 : 300 mm 
Fig. 9 Effects of cohesive force and measurement height on 
the average frequency of pressure signal. Fluidization gas 
velocity: 2.5Umf, static bed height: 300 mm 


3 结 论 
本 文 结合 X 射线 层 析 成 像 和 压力 脉动 分 析 技术 
研究 黏 性 B 类 颗粒 在 三 维 鼓 泡 流 化 床 中 的 流 化 特性 ， 
获取 了 不 同 黏 性 力 下 的 气泡 形态 演变 ， 流 化 失效 方 
式 以 及 失效 前 后 流 化 特征 参数 的 变化 规律 ， 得 到 主 
要 结论 如 下 : 
1) 当 吾 性 力 较 大 时 , 流 化 以 节 涌 的 方式 失效 , 而 
且 由 于 料 检 的 易 破碎 性 ， 床 内 交 蔡 产生 节 涌 和 正常 
流 化 ， 随 着 茜 性 力 的 增加 ， 流 化 过 程 中 节 涌 持续 时 
长 所 占 比 例 逐 渐 增 加 ; 
2) 黏 性 力 的 存在 促进 了 三 维 床 中 气泡 的 合并 ， 
导致 气泡 频率 逐渐 降低 ， 气 泡 尺 寸 逐 渐 增 加 ， 由 
于 低速 气 栓 在 气泡 中 的 占 比 与 备 性 力 呈 正比 关系 ， 
因此 床 内 的 平均 气泡 上 升 速度 随 弗 性 力 的 增加 逐渐 
降低 ; 
3) 伴随 着 节 涌 的 产生 ， 压 力 信号 的 功率 谱 密 度 
呈 双 峰 分 布 ， 低 频 峰 对 应 节 涌 ， 高 频 峰 对 应 正常 流 
化 . 节 涌 在 压力 信号 能 量 分 布 中 所 占 比例 随 黏 性 力 
的 增加 而 增加 , 导致 平均 频率 在 床 温 超过 35°C 后 迅 
速 下 降 ， 
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